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ABSTRACT 

We present a method of determining lower limits on the masses of pre-main-sequence 
(PMS) stars and so constraining the PMS evolutionary tracks. This method uses 
the red-shifted absorption feature observed in some emission-line profiles of T Tauri 
stars indicative of infall. The maximum velocity of the accreting material measures 
the potential energy at the stellar surface, which, combined with an observational 
determination of the stellar radius, yields the stellar mass. This estimate is a lower 
limit owing to uncertainties in the geometry and projection effects. Using available 
data, we show that the computed lower limits can be larger than the masses derived 
from PMS evolutionary tracks for M < 0.5M Q . Our analysis also supports the notion 
that accretion streams do not impact near the stellar poles but probably hit the stellar 
surface at moderate latitudes. 
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1 INTRODUCTION 

A complete understanding of the star formation process re- 
quires the ability to predict how the properties of young 
stars depend on the initial conditions of star formation. De- 
spite significant observational effort in the last two decades, 
a key experimental problem remains. That is reliable de- 
termination of the masses and ages of pre-main sequence 
(PMS) stars. Stellar masses are needed to investigate how 
the initial mass function (IMF) varies from region to region 
as a function of initial conditions (i.e. cloud parameters). 
With reliable stellar ages in hand, we can begin to address 
quantitatively important questions of early stellar evolution. 
Presently, a common method used to determine masses and 
ages of young stars is to place them in an H-R diagram and 
compare their positions with theoretical evolutionary tracks 
(eg Hillenbrand 1997). Considerable theoretical work has 
been done recently that has advanced our understanding of 
PMS evolution (D'Antona & Mazzitelli 1994 [DM94]; Swen- 
son et al. 1994). Yet there remain significant differences 
in the tracks owing to alternative treatments of convection, 
opacities and the zero-point of the calculated ages. An ad- 
ditional complication is on-going accretion during the PMS, 
which can alter the path of a star during its early evolu- 
tion (Hartman, Cassen & Kenyon 1997; Seis, Forestini & 
Bertout 1997). Without adequate observational constraints, 



it is nearly impossible to determine what the correct treat- 
ment should be. 

The discrepancies in the PMS tracks are largest for the 
lowest-mass stars (Hartigan, Strom, & Strom 1994). It is 
precisely this end of the stellar mass distribution which is the 
most uncertain both amongst field stars as well as in star- 
forming regions (eg Kroupa, Tout, & Gilmore 1993; Luh- 
mann & Rieke 1998). Bona fide brown dwarf objects have 
recently been discovered as companions to field stars (Naka- 
jima et al. 1995) and as free-floating members of young 
clusters (Rebolo, Zapatero-Osorio & Martin 1995). How- 
ever the frequency of such objects remains unknown. Young 
clusters found in regions of star-formation provide one of the 
best opportunities to determine the relative contribution of 
low-mass stars and brown dwarfs to the cluster IMF because 
substellar objects cool as they age, becoming extremely dif- 
ficult to detect when old. Uncovering local maxima or min- 
ima in the IMF near the stellar/substellar boundary would 
point to the existance of characteristic masses in the forma- 
tion process (Adams & Fatuzzo 1996). Because of the im- 
portance of the low-mass end of the IMF in star-formation 
theory, it is imperative to establish reliable PMS evolution- 
ary tracks and thus PMS masses. To do this, we need in- 
dependent determinations of PMS masses to constrain the 
tracks. 
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The best techniques available for directly determining 
stellar masses involve the study of binary star systems. 
Eclipsing binary systems provide the best estimates of stel- 
lar masses as well as radii (Popper 1980). Unfortunately, 
only a few such systems are known that are sufficiently 
young to provide a test of the PMS evolutionary tracks (eg 
Casey et al. 1998) and even these do not probe the low- 
est masses. Astrometric determination of orbital parame- 
ters for visual PMS binary star systems is another promis- 
ing technique (Ghez et al. 1995), but useful constraints are 
still several years away. Dynamical mass ratios can be de- 
rived for double-lined spectroscopic binaries and compared 
with theoretical mass ratios in order to test the PMS tracks 
(Lee 1992). However most of the known double-lined PMS 
SBs have primary stars with masses greater than solar and 
mass ratios near one. We therefore need other constraints, 
even if they can provide only lower or upper limits on the 
stellar mass. Such limits can be obtained from the dynam- 
ical information available from the material accreting on to 
T Tauri and other young stars. 

Accretion on to young stars occurs through a viscous 
circumstellar disc, where angular momentum transport out- 
wards permits the inward mass transport (Pringle 1981; 
Lin & Pringle 1990). There is a growing body of observa- 
tional evidence indicating that the accretion of material on 
to young stars occurs via magnetospheric accretion columns 
which extend several stellar radii from the PMS photosphere 
to the inner edge of a circumstellar disc (Edwards 1997). 
Material at this point falls freely along the field lines on 
to the stellar surface. Evidence for this picture comes pre- 
dominantly from the observed redshifted absorption compo- 
nent (several hundred kms -1 ) seen in the line profiles of 
the higher Balmer, Pa/3, Br7, Hel and NaD emission lines 
(Walker 1972, Appenzeller, Reitermann & Stahl 1988; Ed- 
wards et al. 1994; Hartmann, Hewett & Calvet 1994; Folha, 
Emerson & Calvet 1997, Folha & Emerson 1998). Models 
of the IR colours of T Tauri stars imply that the circumstel- 
lar discs have inner holes comparable (but slightly interior) 
to the magnetospheric radius which is typically estimated 
at 6-8 R* (Kenyon, Yi & Hartmann 1996; Meyer, Calvet 
& Hillenbrand 1997; Armitage, Clarke & Tout 1998). The 
magnetospheric accretion paradigm also provides an attrac- 
tive explanation of the relatively slow spin rates of T Tauri 
stars (Bouvier et al. 1993; Edwards et al. 1993). The stellar 
magnetic field interacts with parts of the disc that are spin- 
ning slower than the star (at radii beyond corotation) and 
the resulting magnetic braking of those stars with discs (the 
actively accreting CTTS) can explain why they have longer 
periods than those (the non-accreting WTTS) stars without 
circumstellar discs (Cameron & Campbell 1993; Armitage 
& Clarke 1996). 

In this paper, we investigate how the observed red- 
shifted absorption components, known as inverse P-Cygni 
(IPC) profiles, can be used as a measurement of the depth 
of the potential well of the star and hence, combined with 
an observational determination of the stellar radius, yield an 
estimate of a lower limit on the stellar mass. Section 2 de- 
scribes the method while Section 3 provides examples, using 
available data, of the mass estimates from this method and 
compares them with published estimates based on the PMS 
tracks. Section 4 shows how the lower limits on the stel- 
lar mass could in principle constrain the accretion stream 



geometry. A discussion and summary is presented in § 5. 



2 CONSTRAINING PMS MASSES 

The IPC profiles in the emission line spectra of classical 
T Tauri (CTTS) stars offer a direct measurement of the 
depth of the potential well of the star. The large infall ve- 
locities seen in the absorption profiles are generally assumed 
to arise from material in free-fall from several stellar radii. 
This near radial infall is believed to be caused by the dis- 
ruption of the accretion disc by the stellar magnetic field 
(Konigl 1991, Edwards et al. 1994). The velocity of the 
stream, as it impacts the stellar surface then directly mea- 
sures the difference in potential energy from the radius at 
which the disc is disrupted to the stellar surface: 

1 2 _ GM± GM± 

„ ^stream — D D , (,1) 

1 ti-k Khole 

where V s tream is the velocity of the accretion stream at im- 
pact, Mi, is the stellar mass, -R* is the stellar radius, Rhole is 
the magnetospheric radius at which the disc is disrupted and 
G is the gravitational constant. Knowing these quantities, 
we can evaluate the stellar mass as 

M* = ( 1 - -ft-) \ (2) 

ZLr \ -Kholc / 

The rotational energy of the disc matter at Rhole is as- 
sumed to be dissipated as the disc matter couples on to 
the star's magnetic field lines, and is therefore neglected in 
equation (1). This energy is at most half the potential en- 
ergy at -Rholc an d thus its contribution to equation (1) is 
equivalent to the disc being disrupted at 2 -Rholc- Mass esti- 
mates incorporating the rotational energy will therefore lie 
in between the two limits we calculate below; i) that the 
matter falls in from infinity; or ii) from -Rholc- 

Although the value of R*/.Rhoic is uncertain, it can be 
estimated by comparing the infrared excess emission ex- 
pected from a circumstellar disc with and without an in- 
ner hole. Kenyon, Yi & Hartmann (1996) estimate that 
Rhole / R* ~ 4 (with acceptable values between 3 and 5). 
Similar results were obtained by Meyer et al. (1997) who 
estimate that 2 < Rhole/R* < 6 for a sample of T Tauri 
stars located in the Taurus dark cloud. In practice, because 
we are looking for a lower limit on the stellar mass, and as 
R*/Rhoic is relatively small, we calculate the stellar mass 
using equation (2) assuming either that the material falls in 
from infinity or from five stellar radii (-Rholc ~ 5R*). We 
note that the case of infall from infinity provides a strict 
lower limit to the mass resulting in values 0.8 times those 
obtained in the 5-R* calculation. 

We see from equation (2) that we can make an estimate 
of the stellar mass from just the stellar radius and the impact 
velocity of the accretion stream. The stellar radius can be 
calculated directly from the stellar luminosity and T e g. It 
is thus independent of the PMS tracks, but does depend 
on observers' ability to correct for reddening, establish the 
spectral type (i.e. photospheric temperature) , and separate 
stellar from accretion luminosity. The velocity itself can 
be estimated from the red-shifted absorption profile. There 
is an added complication that the accretion stream is not 
necessarily parallel to our line of sight so that there is a 
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projection effect in the velocity. The observed IPC velocity 
relates to the impact velocity as 

Vipc = Stream X COS0, (3) 

where 8 measures the angle between our line of sight and 
the direction of the stream on impact. This implies that 
the observed Vipc will be a lower limit on the true impact 
velocity and again gives us a lower limit on the stellar mass. 
Furthermore, as the accretion streams most probably rotate 
with the star, they will be obscured by the star over part of 
the rotation period (eg Smith et al. 1997) and the observed 
V; pc wm be variable (Edwards et al. 1994). 

In the above we have not considered the exact mecha- 
nism by which the disc matter couples to the magnetic field 
and the accompanying torques. This simplification neglects 
any azimuthal motions of the infalling matter as it follows 
the field lines. Fortunately, any such deviations also ensure 
that the mass estimates are lower limits. In fact, all of the 
unknowns involved in measuring the star's gravitational po- 
tential well, the distance from where the material falls, the 
projection of the stream along our line of sight, that the 
stream impacts radially, and even the assumption that the 
stream's maximum velocity is on the stellar surface and not 
at some distance from the star, make our estimate a lower 
limit. Therefore we should, in general, underestimate the 
gravitational potential well of a star. Given an accurate 
determination of the stellar radius from its luminosity and 
temperature, the mass determined from this method is a 
lower limit on the true stellar mass. 



3 COMPARISON OF IPC MASS LIMITS AND 
MASSES FROM PMS TRACKS 

In order to ascertain how useful this method is for constrain- 
ing PMS masses, we compare masses derived from PMS 
tracks (M trac k s ) with the lower limits on the stellar masses 
calculated with equation (2) (M; pc ). Table 1 shows this 
comparison for 13 CTTS stars chosen from the stars stud- 
ied by Hartigan et al. (1995). This sample includes all stars 
with observed IPC profiles (from Edwards et al. 1994, Folha 
& Emerson 1998 and in one case Appenzeller et al. 1988) 
that have also been surveyed and found not to have any 
companions within separations 0.1 — 2.0" (Mathieu 1994 
and references therein; Simon & Prato 1995). Companions 
in this separation range are unresolved in the Hartigan et 
al. (1995) study and hence would introduce errors into the 
determination of the stellar luminosity and thus the radii 
and PMS track masses (Chez, White & Simon 1997). This 
highlights another advantage of the technique; it can be ap- 
plied to single stars and not just binary systems in con- 
trast to the methods for determining directly stellar masses 
discussed above. The radii and masses (from D'Antona & 
Mazzitelli 1994 tracks adopting Alexander opacities and CM 
convection) are taken from Hartigan et al. 1995. These au- 
thors have attempted to deredden the stars and remove the 
effects of accretion luminosity before calculating the stellar 
parameters. 

The velocities V5 pc are estimates based on the IPC pro- 
files seen in the Balmer, NaD and Hel lines (Edwards et 
al. 1994; Appenzeller et al. 1988), and in the Paschen/3 and 



Brackett7 lines (Folha & Emerson 1998). The use of differ- 
ent lines to measure the impact velocities is non-ideal as it 
may introduce uncertainties due to radiative transfer effects, 
but is necessary because of the paucity of available IPC pro- 
files. For each star, two estimates of the velocity were made, 
the first from the deepest part of the absorption feature and 
the second from the maximum velocity of the absorption 
(taken to be where the absorption profile meets the contin- 
uum). These estimates are most probably lower limits on the 
velocity of the infalling material when it impacts the stellar 
surface. Modelling of the IPC profile with radiative transfer 
(Muzerolle, Calvet & Hartmann 1998) shows that the im- 
pact velocity typically corresponds to the red-most extreme 
part of the absorption component. The minimum in the IPC 
profiles in Muzerolle et al. (1998) generally corresponds to 
a Vip C that is 10 to 20 per cent lower than the character- 
istic velocity V s tream- We prefer to use both estimates of 
Vipc (quoted in Table 1) as limits of the true characteristic 
velocity to avoid any uncertainties in the details of the line 
formation. For each star, the maximum velocity observed (of 
each type) was used. In general, IPC profiles are variable (eg 
Edwards et al. 1994; Smith et al. 1997) with higher velocity 
IPC profiles presumably arising when the accretion stream 
is most closely aligned to our line of sight. These veloci- 
ties are therefore lower limits on the actual impact velocity 
of the accretion stream. In this context it is worth noting 
the two different estimates of the IPC velocity profile of DR 
Tau (from Edwards et al. 1994 and Appenzeller et al. 1988) 
in Table 1. The IPC profile can even disappear completely 
when the stream is behind the star. This implies that the 
magnetospheric accretion geometry is complex and is not 
an axisymmetric ring aligned with the rotation axis. The 
dynamical mass estimates (lower limits) based on the IPC 
velocity measurements are given in columns 5 and 6 of Ta- 
ble 1 assuming that the stream velocity, Vstream, is given by 
the position of the minimum of the absorption feature, and 
that the matter falls in from infinity (-R*/-Rhole = 0) or from 
five stellar radii (R*/Rholc = 0-2). The corresponding cases 
for the maximum velocity seen in the IPC profile are given in 
columns 8 and 9. The photometric (rotation periods, v sini, 
and inclinations of the stars, derived by Bouvier 1995, are 
given in column 10. 

It is difficult to quantify the uncertainties in our mass 
determinations presented in Table 1. First we consider the 
errors in the stellar radii derived from the observational 
determination of the stellar luminosity and effective tem- 
perature. Propagating errors in the photometry, redden- 
ing, distance modulus, bolometric corrections (i.e. luminos- 
ity), and spectral types (i.e. photospheric temperatures), 
we estimate that the radii are accurate to within 20-30 per 
cent (in agreement with the analysis of Kenyon & Hart- 
mann 1995). We can estimate the systematic uncertainties 
in the radii by comparing the derived values from differ- 
ent observational determinations of the stellar parameters 
(Hartigan et al. 1994; Gullbring et al. 1998). Both stud- 
ies rely on independent spectra and use different corrections 
for the photometry. The stellar radii thus derived indepen- 
dently agree remarkably well, with differences at the 10-15 
per cent level. 

Uncertainties in our estimates of Vi pc are approximately 
10 per cent. Furthermore, we neglect radiative transfer ef- 
fects and the errors associated with estimating the maximum 
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Table 1. Comparison of T Tauri masses from PMS tracks and lower limits from IPC profiles. The radii (in Rq), and PMS track 
masses (in Mq) are taken from Hartigan et al. 1995 using the tracks from D'Antona & Mazzitclli 1994. The IPC profile velocities (in 
km s _1 ) are from 1 Edwards et al. 1994, from 2 Folha & Emerson (1998) and from 3 Appenzeller et al. 1988. Two values of V;p C are 
used: these are the deepest absorption, V;p C (deep) corresponding to M; pc , and the the maximum velocity in the absorption profile 
Vjp C (max) corresponding to M; pc * (see text). The masses, M; pc /Mq, are calculated on the assumptions that the matter falls in from 
infinity, Mjp C (oo), and from 5-R*, M;p C (5i?*). The last three columns list the rotational periods in days, vsvcii in km s , and the 
inclination angles (from Bouvier et al. 1995), adapted to the stellar radii from Hartigan et al. 1995). An inclination angles of 90deg is 
assigned to CM Aur where the above comparison yields sini > 1 



Star 


Radius 


M tracks 


line Vj 


pc (dcep) 


M ipc (oo) 


M ipc (5fe) 


y ipc (max) 


Mip C (oo)* 


M ipc (5i?*)* 


Prot 


vsini 


i 


AA Tau 


1.8 


0.38 


H/3 


220 1 


0.23 


0.28 


280 1 


0.34 


0.41 


8.2 


11.4 


81 


BP Tau 


1.9 


0.45 


H<5 


290 1 


0.42 


0.51 


400 1 


0.8 


0.96 


7.6 


7.8 


38 


CW Tau 


2.2 


1.03 


Br7 


170 2 


0.17 


0.20 


260 2 


0.40 


0.48 








DK Tau 


2.7 


0.38 


H 7 


280 1 


0.56 


0.67 


315 1 


0.70 


0.84 


8.4 


11.4 


43 


DL Tau 


1.9 


0.37 


Na D 


240 1 


0.29 


0.34 


315 1 


0.50 


0.60 








DN Tau 


2.2 


0.42 


H/3 


160 1 


0.15 


0.18 


205 1 


0.25 


0.30 


6.0 


8.1 


27 


DO Tau 


2.4 


0.31 


Pa/3 


175 2 


0.19 


0.23 


235 2 


0.34 


0.41 








DR. Tau 


2.7 


0.38 


Hel 


225 1 


0.36 


0.43 


315 1 


0.71 


0.85 














H<5 


330 3 


0.77 


0.93 














DS Tau 


1.6 


1.28 


H<5 


265 1 


0.29 


0.35 


370 1 


0.56 


0.68 








FM Tau 


1.6 


0.15 


Pa/3 


120 2 


0.06 


0.07 


210 2 


0.18 


0.21 








CI Tau 


2.5 


0.30 


Pa/3 


225 2 


0.33 


0.40 


350 2 


0.81 


0.97 


7.2 


11.2 


40 


GK Tau 


2.2 


0.41 


H/3 


250 1 


0.36 


0.43 


280 1 


0.46 


0.55 


4.65 


18.7 


52 


CM Aur 


1.6 


0.52 


Pa/3 


250 2 


0.26 


0.31 


315 2 


0.41 


0.49 


12.0 


12.4 


90 



absorption velocity in a noisy spectra. Both of these should 
tend to reduce the observed velocity from the true impact 
velocity. Thus, as we are deriving lower limits, uncertain- 
ties in our estimates of V5 pc from the data should be smaller 
than the difference between our two limits. Combining the 
errors in the radii with the 10 per cent error in estimating 
the velocities, this corresponds to a maximum uncertainty 
in our lower limits of ±45 per cent. 

Of course the errors in deriving the stellar luminosity 
and temperature also result in uncertainties in the mass es- 
timates from the PMS tracks. Because the mass tracks arc 
generally vertical during the Hyashi contraction phase for 
late-type stars, errors in spectral type and conversion to ef- 
fective temperature dominate the errors in estimating stel- 
lar masses from PMS evolutionary models. Typical errors 
of ±1 spectral subclass translate into errors of ±0.02 dex 
in log(T cf f), resulting in relative errors of ±O.1M0 in stel- 
lar mass estimates for young (< 10 Myr) late- type (K-M) 
stars for the DM94 tracks (Hillenbrand 1997). Apart from 
problems with the tracks themselves (which we attempt to 
probe with this technique), different approaches for placing 
stars in the H-R diagram taken by different groups could 
be an additional source of systematic error. The masses re- 
ported by Gullbring et al. (1998) (derived from the same 
DM94 tracks) generally agree with those reported by Harti- 
gan et al. (1995) except in those cases where a significantly 
different spectral type was adopted. Gullbring et al. 1997 
adopted a later spectral type for DS Tau than HEG (K5 vs. 
K2) and an earlier spectral type for GI Tau (K6 vs. M0). 

The ratio of the dynamical mass estimate, Mj pc , to the 
track mass, M^a^s, is plotted in Figure 1 against the PMS 
track mass and the radius. The ratio is plotted for both 
mass estimates based on the two characteristic velocities and 
on the hole size in the disc (matter falling in from infinity 
(filled) and from 5 stellar radii (open)). As the dynami- 



cal mass estimates are lower limits on the stellar masses, 
the ratio plotted should be less than one for all stars. We 
see that this is not the case for some of the stars when we 
use the lower characteristic velocity (upper panels) and with 
most stars if we use the maximum velocity seen in the IPC 
profile (lower panels). The three stars for which this dis- 
crepancy is most striking (DK Tau, DR Tau, and GI Tau) 
are all unusual and have larger radii than the other stars. 
They should be treated with some caution. DK Tau is a 
relatively wide binary (2.5" separation) which should, in 
principle, be uncontaminated by it's companion. Contri- 
butions from a companion can lead to a larger luminosity 
(radius) and hence mass using the IPC method while not 
significantly effecting the mass derived from the tracks. GI 
Tau has M; pc < M tracks according to the stellar parame- 
ters of Gullbring et al. (1998) while the high veiling of DR 
Tau makes it's stellar parameters somewhat uncertain. The 
other stars commonly have mass limits that are similar to 
or lower than the track mass when the lower estimates for 
Vj pc are used. If we use the maximum velocity seen in ab- 
sorption, then almost all the derived dynamical mass limits 
have M ipc > M tracks . 

Furthermore, some of the stars should be seen with rel- 
atively large projection angles, so we would expect them 
to have M; pc < 0.7M trac k s (for a projection of 45 degrees). 
Thus, from the fact that most of the stars with masses 
M* < O.5M0 have M; pc > 0.7M tracks , even when the veloc- 
ity is estimated from the deepest part of the absorption, we 
deduce that the DM94 tracks potentially underestimate the 
true stellar masses. This problem is aggravated if the disc is 
truncated at smaller radii, resulting in higher estimates of 
M; pc . Therefore, we can conclude that the dynamical mass 
estimates provide a useful lower limit on the stellar mass 
and thus a constraint on the PMS evolutionary tracks. 
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Figure 1. The ratio of the dynamical mass estimate, M; pc , to the mass derived from the PMS tracks, M trac k s , is plotted against 
^tracks / M@ (left panels) and against the stellar radii in Rq (right panels). The top panels (circles) use the deepest part of the 
absorption profile as being the characteristic V;p C , while the bottom panels (squares)assume that V;p C is the maximum velocity in the 
absorption profile. The filled symbols assume the matter falls in from infinity while the open symbols assume that the mass falls in from 
five stellar radii. As Mjp C is a strict lower limit, the ratio of M;p C /M trac k s should always be less than 1 (dashed line) if the PMS track 
masses are correct. The dotted line represents the expected ratio for a projection of 45 degrees. 



Other PMS tracks give different results. For example, 
the tracks of Swenson et al. (1994) give systematically higher 
masses (see Hillenbrand 1997 for comparison) which would 
compare more favourably with our upper-limits. However 
this technique can in principle only reject a set of tracks 
for which the predicted masses were too low. As long as 
the track masses are comfortably above the lower limits 
set by Mj pc we cannot prefer one set of tracks over an- 
other. Nonetheless, a direct comparison of different PMS 
track masses (including the effects of accretion) with mass 
limits derived using the method described here for a statisti- 
cally significant sample of PMS stars will assist in discarding 
unviable PMS evolutionary tracks. 

4 INCLINATION AND ACCRETION STREAM 
GEOMETRY 

The dynamical mass estimate of equation (2) can also tell 
us something about the geometry of the magnetospheric ac- 



cretion streams. The mass estimate uses the observed infall 
velocity along our line of sight and so can constrain the ori- 
entation of the stream. Once reliable masses are available 
from the PMS tracks we can combine the mass estimate 
from the IPC velocity with the knowledge of the star's incli- 
nation to our line of sight and constrain where on the star 
the accretion stream impacts. Figure 2 plots the ratio of the 
dynamical mass estimate, M; pc (using the lower velocity cor- 
responding to the deepest part of the absorption profile), to 
the mass derived from the PMS tracks, M trac k s , against the 
inclination angle of the star for those stars with both mea- 
sured rotation periods and wsini determinations (from Bou- 
vier et al. 1995,). The star GI Tau is excluded from the plot 
because of the uncertainty in its radius (between Hartigan et 
al. [1995] and Gullbring et al. [1998]) and hence the inclina- 
tion. The inclination angles are also listed in Table 1. GM 
Aur is assigned an inclination of 90deg because the apparent 
sini > 1. From Figure 2, we see that stars oriented pole-on 
(0 deg inclination) and equator-on (90 dcg inclination) have 
dynamical mass estimates significantly lower than those de- 
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Figure 2. The ratio of the dynamical mass estimate, M; pc , to 
the mass derived from the PMS tracks, Mj. rac k s , is plotted 
against the star's inclination to the line of sight (pole-on 
corresponds to degrees). The dynamical mass, M; pc , is 
calculated from the lower characteristic velocity representing the 
deepest absorption and the assumption that matter falls in from 
infinity (filled circles) or from five stellar radii (open circles). 
The smaller symbols indicate an uncertainty in the inclination 
of less than 10 degrees while the larger symbols represent an 
uncertainty in the inclination of less than 5 degrees. The highest 
point corresponds to DK Tau (a 2.5" binary). GI Tau has not 
been plotted because of the large difference in radii between 
Hartigan ct al. (1995) and Gullbring et al. (1998), which leads 
to a large uncertainty in the inclination. 

rived from the PMS tracks. In contrast, stars with moderate 
inclination angles (30 to 60 deg) have dynamical mass esti- 
mates approximately equal to or greater than those derived 
from PMS tracks. 

Although the sample size is small, Figure 2 hints that 
the accretion streams may impact preferentially at moderate 
latitudes. This argues against scenarios of magnetospheric 
accretion where the accretion streams are assumed to impact 
on, or near, the poles. In contrast, Muzerolle et al. (1998) 
adopt a dipolar geometry where the accretion stream im- 
pacts at moderate latitudes but that implies small magne- 
tospheric radii such that Rholc ~ 2.R*. Kenyon et al. 1994 
derive a similar model based on extensive multi-color ob- 
servations of the rotating "hot spot" observed on DR Tau. 
A pure aligned dipolar field implies a relation between the 
inclination at which the accretion stream impacts (<f>, mea- 
sured from the pole) and the ratio of Rh i c /R* given by 
(Hartmann et al. 1994) 

-Rholo _ - -2 , 



(4) 



As the observed Vi pc for the accretion stream should have a 
maximum where the star's inclination i is equal to this angle 
4>, then we would expect the inclination angles of 40 to 60 
deg (as seen in Figure 2) to correspond to disc truncation 
radii of Rholc = 1-33 to 2.4.R*, much smaller than the 5 R* 
used here. These smaller values would in turn imply signifi- 
cantly larger values of Mi pc . This may imply that accretion 
geometry is complex, and not a pure, aligned dipole. Using 



the method described here, with a large enough sample, it 
should be possible to constrain the geometry of the accretion 
region. 

5 DISCUSSION AND CONCLUSIONS 

Determining accurate masses for large ensembles of pre- 
main-sequence (PMS) objects is presently impossible owing 
to our inability to adequately constrain the PMS evolution- 
ary tracks. This limitation is a major obstacle in determin- 
ing the low-mass end (0.05 — 0.5Mq) of the IMF in regions 
of recent star formation (Hillenbrand 1997). Here, we have 
shown that lower limits on the stellar masses can be deter- 
mined by measuring infall velocities on to these objects. 

In summary, the presence of red-shifted absorption pro- 
files in the emission line spectrum of PMS stars indicates 
near radial accretion. In the magnetospheric accretion 
model, the stellar magnetic field disrupts the circumstel- 
lar disc at several stellar radii from where the matter falls 
freely in along the field lines on to the stellar surface. The 
velocity of the infalling matter (from the IPC profile) mea- 
sures the potential energy at the stellar surface which, when 
combined with a determination of the stellar radius, yields 
an estimate of the stellar mass. 

For a collection of 13 stars (Hartigan et al. 1995), 11 
have lower limits which are comparable to or greater than 
the track masses and thus can act as significant constraints 
on the stellar mass. Using the deepest part of the absorp- 
tion profile as the impact velocity, and assuming the disc 
is truncated at five stellar radii (-Rholc ~ 5R*), 5 of the 13 
stars have M; pc > M trac i; S . If the disc is truncated at in- 
finity, three of these five stars still have M; pc > M tracks . 
If we model the impact velocity as the maximum velocity 
seen in absorption, then 10 of the 13 stars have lower lim- 
its M ipc > M tracks , 8 of which have M ipc > M tracks if the 
disc is truncated at infinity instead of at 5R*. Truncation 
radii smaller than 57?* result in larger dynamical mass lim- 
its such that if -Rholc ~ 2R* , none of the lower limits would 
be compatible with the PMS tracks. From this we can con- 
clude that i) the method described here is a useful method 
for estimating PMS masses and thus constraining the PMS 
evolutionary tracks, and ii) that there are potentially signif- 
icant problems in the DM94 tracks. 

In addition, this method for setting lower limits on the 
stellar mass can also help us constrain the geometry of the 
magnetospheric accretion. Using information on the star's 
inclination to our line of sight, we find that stars seen at 
moderate inclination appear to suffer less from projection 
effects. This implies that the accretion stream probably im- 
pacts on the star at moderate latitudes. 

Although this technique appears promising, additional 
work is required before quantitative tests of the available 
PMS tracks can be undertaken. Specifically, we require; i) 
more extensive collections of IPC profiles, ii) improved es- 
timates of stellar radii based on careful photometric and 
spectroscopic observations and iii) additional modelling ef- 
forts to determine the appropriate velocity to adopt from 
the IPC line profile, in order for this method to become a 
powerful tool in the constraint of the masses of PMS stars. 
With these data, we can then begin to examine the PMS 
evolutionary tracks and study the geometry of the accretion 
streams. 
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